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1. Introduction 

It is well known that theories with general covariance are constrained systems [1- 
4]. In other words, their equations of motion in Lagrangian formalism lead to 
acceleration-free relations. On the other hand, constructing the Hamiltonian for- 
mulation for such theories needs care in order to consider the constraints. The 
primary constraints emerge in the phase space whenever the momenta are not inde- 
pendent functions of velocities. The secondary constraints come out as the result of 
the consistency of primary constraints. 

The most difficulty in Hamiltonian treatment of general covariant theories is that 
the action depends on the second derivatives of the metric, as well. In this situation 
a wellbehaved Hamiltonian system is not recognized, or at least is not agreed upon, 
even when a system is not constrained. However, for Einstein-Hilbert gravity [5] or, 
for instance, Hofava gravity [6] one may use one may use the Arnowitt-Deser-Misner 
variables [7] whose advantage is that the Lagrangian does not contain accelerations 
when written in terms of these variables. This may not happen for an arbitrary 
general covariant Lagrangian. 

The other possibility is using the so-called Palatini formalism in which the 
Christoffel symbols are considered as independent variable. For Einstein-Hilbert 
gravity, this approach does work well. The reason is the relation between the 
Christoffel symbols and the derivatives of the metric results naturally from the equa- 
tion of motion of Christoffel symbols. It was shown recently [8] that for a gravitation 
theory of the Lovelock-type, the Palatini formalism is fine. However, for an arbitrary 
model the equations of motion give no guaranty about the relations of Christoffel 
symbols and derivatives of the metric. Therefore, one needs to add them to the 
Lagrangian using Lagrange multipliers, which should be considered as additional 
variables in the Lagrangian formalism. 



- 1 - 



The three-dimensional gravity has attracted intensive interest in recent years. 
One reason is that it is possible to construct nontrivial renormalizable models in 
three dimensions. Among so many attractive features, investigating the Hamiltonian 
structure of the models is noticeable. The topological massive gravity [9], which is 
generally covariant on a closed manifold, is one of the most important ones. The 
Hamiltonian structure of TMG is discussed in some papers [10]. 

Recently Begshoeff, Hohm and Townsend proposed a model [11] for three dimen- 
sional massive gravity (the so called new massive gravity) which preserves parity and 
possesses general covariance on an arbitrary manifold. Linearization around the flat 
metric of this model leads to Pauli-Fierz action describing massive graviton. Then 
Deser [12] showed that this model is finite and ghost-free. Oda and Nakasone [13] 
showed afterward that the model is unitary and renormalizable. Clement also gave 
some black hole solutions of the model [14]. The Lagrangian of new massive gravity 
(NMG), at one hand, includes accelerations (i.e. second order derivatives of the met- 
ric), which make it necessary to use Christoffel symbols (or combinations of them) as 
auxiliary fields. On the other hand, the model is not of the Lovelock-type. These pe- 
culiarities make the Hamiltonian treatment and constraint structure of NMG much 
more difficult. Moreover, the existence of quadratic terms with respect to makes 
it difficult to use the Arnowitt-Deser-Misner variables. However, in spite of compli- 
cated calculations, the Hamiltonian treatment of the theory can be followed carefully. 
This is what we have done in this paper. 

Our main task in this work is counting the physical degrees of freedom of this 
model. Since the Einstein gravity in three dimensions has zero degrees of freedom, one 
expects roughly that the NMG action, which contains two more derivatives, should 
have 2 degrees of freedom; in the same way as TMG with one more derivative than 
Einstein- Hilbert action possesses 1 degree of freedom. This result is in agreement 
with the number of massive gravitons in the linearized model. However, from a 
theoretical point of view, it is better to check the validity of such rough arguments 
by a careful Hamiltonian analysis. 

In Section 2 we introduce the model in the Hamiltonian formalism and find 
the primary constraints. In Section 3 we follow the consistency conditions of the 
constraints and find the secondary constraints of system. Section 4 is devoted to our 
conclusions. 

In our work we use Greek indices for space-time components and Latin indices 
for space components. 
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2. Lagrangian and Hamiltonian 



The action of NMG is given as: 

1 ,„ „„„ 3 



R-2A+ -^(R^Rf - r B 2) 



rrr 



(2.1) 

where g is the metric determinant, R pu is the Ricci tensor and R is the Ricci scalar. 
We assume that 3d space-time is torsion-free and the Christoffel symbols are 
symmetric with respect to \i and v. This allows us to introduce new variables 
via 



= ££/ - + S u^°la), (2-2) 

as in reference [15]. The Ricci tensor in terms of £ variables contains derivatives in 
the form of total 3-divergence as 

-^V = £fjtv,\ ~ £fj,a£v\ + 2^mo-Ca- (2-3) 

This is, in fact the advantage of using £'s in comparison with T's. In this way, ^ u 
are the only variables whose velocities are present in the Lagrangian. As is well 
known, in 3 dimensions an action containing higher order derivatives can not be of 
Lovelock-type. So the Palatini approach can not be used without imposing explicitly 
the relation between the metric g^ v and auxiliary variables . Using the Eq. (|2.2| ) 
and the definition of Christoffel symbols as 

= + 9vw - Qwp), (2-4) 



we have 



= 9x P (C(s - + - \{gxp,a + 9xa,p - g a 0,x) = 0. (2.5) 



Expressions ^ a \p should be considered as external Lagrangian constraints which 
should be put by hand in the Lagrangian with Lagrange multipliers. These Lagrange 
multipliers then should be taken into account as new variables in addition to g^ v and 
. In this way, the following action should be considered instead of the original one 

& 

S = J d 3 xC = J d 3 x^g-[R -2A + ^{R^BT - ^R 2 )} 

+ J d 3 xA aX ? ( g Xp (C, - \{5 p a ^ + " \{9xp,a + 9x*,p - g a p,x)) (2.6) 
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where A aX/3 is Lagrange multiplier. Let us enumerate different degrees of freedom 
before considering the details of the dynamics of the system. We have six degrees of 
freedom g^ u and 3 x 6 = 18 auxiliary variables , taking into account the n -H- v 
symmetry in both cases. In other words, by adding the auxiliary variables, number of 
Lagrangian degrees of freedom are multiplied by 4 to avoid higher order derivatives 
(greater than 2) in the equations of motion. We have also introduced 18 more degrees 
of freedom A aX/3 , since Eq. ( |2.5| ) is symmetric with respect to indices a and /? and we 
have A " = A^ Xa . Putting all these points together we have a priory 18 + 18 + 6 = 42 
Lagrangian variables which is equivalent to 84 phase space variables. 
Now we proceed to the Hamiltonian formalism. The canonical momenta conjugate 
to g^ u , and A aX/3 are defined respectively as 

3C 2 

K = = v^or + —G^R a ,), (2.7) 



7T^ 



3C 1 

-(A°^ + A°^-A^), (2. 



dg^u 2 



dL 

PaXB = —■ = 0. (2.9) 

dA aX P 



where G^ avl3 is the generalized metric defined by 



G Hoa,p = g m g v$ _ -g^g^. (2.10) 



From eq.(|2.7|), with 7 = 0, we have 



itW = Im'g^ig-^Uf - (2.11) 



where Q m v$ is the inverse of the generalized metric C" 1 ^ such that Q^avpG^ 
5°5p. Hence, using ( |2.3| ) we find 

& = \ m2 9~ 1/2 + 4m 2 ^ - ? mi + QMx ~ \%Aix- ( 2 - 12 ) 



Using eqs.( p.ll"D and ( 2.12 ) the canonical Hamiltonian density can be derived in the 
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usual way as He = J d 2 xHc-, where 

U c = vr^ V + ICO, + P«xpA aX P - £ 



mV 1/2 ^nW + 2v^(A - 3m 2 ) 



+ -(A^ U + A i ^-A^ u )g i 
2 



TT AU/ 



aA /-cr _j_ £-cr /-A 



(2.13) 



In deriving the Canonical Hamiltonian the following primary constraints resulted 
from Eqs. ( fE7| - [2.9|) are imposed: 



^fW _|_ _|_ _^0i/jti _ ^MO^j 



0. 



n 



0. 



(2.14) 



where the symbol "~" means weak equality, i.e., equality on the constraint surface. 
We recall that primary constraints are identities amongst coordinates and momenta 
which follow directly from the definition of Canonical momenta. The total Hamilto- 
nian reads 



Hn 



SxH T 



fiav ; 



(2.15) 



where U^ u , X 1 and V fJ,au are Lagrange multipliers (in the context of Hamiltonian 
constrained systems) corresponding to the primary constraints ( |2.14j ) respectively. 
The fundamental Poisson brackets of field variables are 



{g iaf (x) t ^(y)} = A^\x-y) i 

{&(x),I^(y)} = 6$A$S&>(*-y), 
{A^(x), P^(y)} = 6$A$S&>(x - y), 



(2.16) 



where A# = + tf; 



3. Constraint dynamics and counting physical degrees of free- 
dom 

The number of primary constraints as well as their corresponding Lagrange multi- 
pliers are as follows: 



- 5 - 



# 4>^ v corresponding to U^ v = 6 

# ^f 1 " corresponding to A' = 2 x 6 = 12 

# ^/jAi/ corresponding to V^ Xu = 3 x 6 = 18 

# total primary constraints = 36 

As in all constrained systems, the primary constraints should be valid in the course 
of time. This means that their Poisson brackets with the total Hamiltonian, which 
is responsible for the dynamics of the system, should vanish. This process is the so 
called "consistency of the constraints." 

It is important to remind the reader that at each step of consistency, two main things 
may happen. If a given constraint has non-vanishing Poisson bracket with some 
primary constraints, the corresponding Lagrange multiplier would be determined in 
terms of phase space variables. This is the case when the related constraint is second 
class l . TThe other possibility is that a new constraint emerges as the consistency of 
the given constraint and corresponding Lagrange multiplier is not determined. These 
constraints at different levels of consistency are called second level, third level, and 
so forth, which altogether are remembered as secondary constraints. The process of 
consistency will continue up to the last level in which either a Lagrange multiplier 
is determined (when we have a chain of second class constraints) or the consistency 
is established identically (when the constraints in the corresponding chain are first 
class). Now, we follow the consistency procedure for our problem. 
Consistency of ^f's causes the following expressions vanish: 

xf = {$f, h t } = -i( W - inj"£ A <$? - \{9xiA^ - g Xa A^5») 

+ 2nJ*&) + p <— ► i/ . (3.i) 

Since <&f have vanishing Poisson brackets with all primary constraints, no term 
containing Lagrange multipliers has appeared in Eq. ( |3.1|) . Therefore, consistency of 
12 primary constraints $f give 12 second level constraints xT ■ consistency of \T 
will be investigated afterward. 
For Q^xv we have 

(3.2) 

The last term above includes Lagrange multipliers U^ v . For fi = i, v = j and A = k 
this term vanishes and the consistency of constraints flikj lead to the following second 

1 Remember that a set of constraints are second class if the matrix of their mutual Poisson 
brackets is non-singular. On the other hand, first class constraints have vanishing Poisson brackets 
with all of the constraints (at least on the constraint surface) . 



{Q, Xu , H T } = g Xp [ ^ - -{HX + CJ P ,) 
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level constraints: 



9kj,i 9ij,k) 



(3.3) 



The constraints 0^ are in fact the same as given in Eq. ( |2.5|) . We should 
investigate the consistency of Qik/s in the next level of consistency. Let come back 
to Eq. ( |3.2| ). The Lagrange multipliers U a p have appeared in the last term due to 
Poisson brackets {O^i/, <j) al3 } with one of indices fx or A or v considered as zero. This 
is a (12 x 6) rectangular matrix of rank 6. So it is possible to divide f^Ai/s with at 
least one zero index into two 6-member sets, as follows: 





= ^001 


B 2 i 


= ^002 




= 2^011 + ^101 


B22 


= 2^022 + ^202 


B\2 


= 2 (^012 + ^021 


k B 12 


= f2oi2 — ^021 



102 





= ^000 


Ci = 


= ^010 




= ^020 


Cn 


= ^011 ~~ 


C22 


= r2o22 — 


, Cl2 


= §(^021 



2^101 
2^202 

— ^012 — 2^102) 



(3.4) 



The constraints of the set f^ 1 -* commute (i.e. has vanishing Poisson brackets) with 
a/ g's and the other set, Vl^ constitute a second class system with Q ^'s, so that the 
6x6 matrix {& 2 \ 0} is nonsingular. We can redefine the Lagrange multipliers 
corresponding to the division of the constraints VL^ Xu into the sets VL tk i , fiW and Q( 2 \ 
such that 

This gives 



V x = V 001 
v {1) - V 

T/2 _ T/002 

v {1) - V 



V, 



(D 



Vg = I(2V W1 + V 0U ) 

V 22 = 2 (2y 202 + T/022) 

yl2 = t/012 + t/021 

T//12 _ T/0 12 lT/102 
V (l) - V ~~ 2 V 



(2) 



' t/0 — T/000 

V(2) - V 

yl _ T/010 

T/2 _ T/020 



Vg = 1{2V 022 - V 202 ) 



(3.6) 



V 12 - T/012 + T/021 _ T/ 



If we consider at this point the consistency of M!/ 's, we get 

{<T,H T } = --Wg^g- 1 ' 2 ^^ - 3n^nD + ^V 1 ^^/^^ 
4 16 

- 2m 2 n(T - y/gsTih - 3m 2 ) + ( Q - + 

- \{A^ V - \A^\i + hv ^ - \v^ v ) + /i 4 — y v. (3.7) 

iLi JLi Zi Zi 
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It can be seen easily that the last term in Eq. ( |3.7|) contains Lagrange multipliers 
V(2) corresponding to the constraints in the set which have nonvanishing Poisson 
brackets with (fi^'s. In this way the set of constraints 
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(3. 



constitute a one-level 12-member family of second class constraints which is con- 
structed of two cross-conjugate chains that determine 12 Lagrange multipliers U a p 
and V(2). In this description, we have used the language of ref. [16] in classifying 
the families of constraints. We just remind here that a family of constraints are 
determined as a set of constraints which are resulted as the consistency of a limited 
subset of primary constraints and make a close algebra of Poisson brackets among 
themselves, and with the canonical Hamiltonian. Consistency of the set Vl^ gives 6 
other constraints of the next level as follows: 

= {D* = i* 022 + * 202 *W = % = { - 2 (3.9) 



where ^ a \p are Lagrangian constraints given in Eq. (|2.5|) . We should continue to 
investigate consistency of the above constraints in the next level. This will make the 
meaning of the classification given in Eq. Q3.9Q more clear. As can be seen, second 
level constraints 0^, ^f^\ and \E^ 2 ) are 12 out of 18 Lagrangian constraints (|2.5). 



The remaining 6 Lagrangian constraints correspond to expressions fl2.5| ) including 
Lagrange multipliers U^ u . In fact the equations of motion of gives = U^. 
Putting this into Eq. (|3.2|) gives the corresponding Lagrangian constraints (|2.5| ) for 
the cases which include time derivatives of the metric. 

Now we proceed to the next level by considering the consistency of xf's, @ifc/ s 
and the sets and fy( 2 \ For 6j fcj 's we find 

Tf? 1 
{©ifcj, Ht} = —^-[gko{Gail3j + GajPi) ~ -(gki(GajfiO + GaO/Sj) 

+ 2^ ki ^ ' 1 ~ ^j^OX + 2^jX^a) + 9kj(£iQ,l — 6^r£oA + o^^Ocr)) 

+ 4m 2 (5f fc0 fi'ij - ^(gkigoj + gkjgoi)) 

+ - + foSftUkx - \{U khl + Uuj - U ijtk ) 

+ gki^lj -1(9^1 + g kj Xu)- (3-io) 



Since U^ v are determined previously, the Eq. (|3.10| ) should be considered as equa- 
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tions to find X l mn . It is easy to check that the matrix of coefficients of A^ n 's , i.e. 
{Qikj, 3*1™™} 5 i s nonsingular and these Lagrange multipliers can be determined com- 
pletely. 

Consistency of constraints in the set gives 



{Dj,H T } = -m 2 g 1/2 (>oo(£o/yv + Qjyau) ~ ^gojGo^H 

o 

— o i.9ooKi + 9ojKi ~ 2<7oiAL), j = 1,2 







(3.11) 



{D' 12 ,H T } 



1 2 -1/2 



9 1 9io(Go,x2u + Gz/iOu)^ 



- T^ioAL + g 20 \\ t - 2g u >Q -1< — >2 



(3.12) 



Since Lagrange multipliers A* fc are already determined, Eqs. (j3.11D -( 3.12 ) are three 
equations for seven unknowns V^, V%\, V/lf, and A* 's. Therefore, we should keep 
these equations in mind and wait for finding four other equations which should be 
solved together with ( |3.11 ) and ( |3.12 ) to find the above unknowns. Anyhow, the 
consistency of Di, D 2 and D' l2 do not go further. Hence, we have a 6 member family 
of second class constraints as 



Si 


B 2 


B 12 




1 


I 


D 1 


D 2 


D' 12 



(3.13) 



- 9 - 



Consistency of xt ' s leads the following expression vanish: 

{XT , H T } = (U^(C X0 - £ A Vo)),* - \(9xoA^ - gXtr A°»%] 



2 g 1 ^((Ga\l3i + Gai/3\)5o ~ -(GaX/30 + GaO/3 A ) $i ) fig ^ fl"^ 







m 

4 



V 



- <U& - \tlX)A xp » + g pa {A apX ^ + A ap ^ 0l 

+ ^{A^Uia - A^t/^) + ^(^^ - ^1/^^) + n o i/.(3.14) 

Let us first consider 6 equations concerning the cases \i — j and u = k in Eq. 
(|3.14 ). It can be seen that the 6x6 matrix of coefficients of Lagrange multipliers 



yikj j g nonsingular. Moreover, the corresponding equations do not include the yet 
undetermined Lagrange multipliers A^ and Aq . They include, however, the Lagrange 
multipliers U^ u , V^) and \\ m which are determined previously. Therefore, the Eq. 
( p.!4|) for the cases considered can be used to determine 6 Lagrange multiplier V* fcj '. 
Nonsingularity of the (6 x 6) sub-matrix {x^ fe ,fi/ mn } has an interesting meaning in 
the terminology of ref. [16] on classifying the constraint families. To this end, the set 
of constraints 

Hi j \ ~' *? 

e ikj ^\xl k , (3-15) 

constitute a family of 24-member, 2-level and cross-conjugate second class system in 
which the consistency of constraints of the first row give the constraints of the second 
row, while the constraints at the end of any chain have nonvanishing Poisson brackets 
with the constraints at the top of the other chain. We can check that {Qikj,®min} 
as well as , xT n } vanish. 

Let us come back to Eq. ( p.!4j) and consider the case (/j, = i, v — 0) or (// = 
0,1/ = i). We have four equations in this case again for seven unknowns V^, V^\, 
V^ 2 and A* 's. These equations are in fact, 4 equations which we were expecting, 
after Eq. fl3.12p . Hence, we have 7 independent equations for 7 unknowns. In this 



way the constraints \i anc l their parents constitute an 8- member, 2-level family 
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of second class constraints shown as 



i 

I (3.16) 

oj 



X 



The only remaining case in Eq. fl3.14p is fj, — v — 0. No term containing Aq 
appears in Eq. (|3.14|) . This corresponds to two constraints Xi° f° r which the term 
including Lagrange multipliers A 00 vanishes. Consistency of xT l ea d to third level 
constraints E 00 as 

£ 00 = -2(11*%), - ((g X0 A°™ - g Xu A^)), - 8m 2 ^n o AO + 2(nf £° + nS°&)& 

+ A 0X0 U iX + g at V°° . (3.17) 

Direct calculation shows that {X 00 ,$° } is a nonsingular matrix. Therefore, consis- 
tency of E 00 determines two Lagrange multipliers A 00 and shows that the constraints 
E 00 as well as their parents in the corresponding chain are second class. In this way 
we have derived a 6-member family of second class constraints gathered in three-level 
chains as 

i 

I 

X? ■ (3.18) 

I 

E 00 



Hence, 12 second class constraints in families ( |3.8| ) determine 12 Lagrange multi- 
pliers U a p and Vt2) at first level of consistency; 38 second class constraints in families 
( EyjD , (|3l6l) and ( gl8l) determine 19 Lagrange multipliers Af , A° J ', V^, y ( 2 1} 

and V^'(i) at second level of consistency; and 6 second class constraints in family 
( |3.18p determine 2 Lagrange multipliers A 00 at third level of consistency. We have a 
total of S=12+38+6=56 second class constraints which determine 12+19+2=33 La- 
grange multipliers. There remain 36-33=3 undetermined Lagrange multipliers which 
are VM, V^, and corresponding to primary constraints Bu,B 22 and B 12 given 

in (& 

Let us recall that consistency of primary constraints Bn,B 22 , and B\ 2 gave us the 
second level constraints Du,D 22 and D\ 2 as in ( |3.9| ). Straightforward calculations 
shows that the Poisson brackets of D n ,D 22 , and D\ 2 with the total Hamiltonian 
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vanishes. Therefore, we collect the following constraints: 



B n 


B22 


B 


1 


I 


I 


D n 


D22 


D 



12 

(3.19) 

'12 

as a 6-member, 2-level and first class family of constraints. So the number of first 
class constraints in the form of family ( |3.19j ) is F = 6. 

For the number of dynamical variables, using the famous formula [17] 

D = N — S — 2F, (3.20) 

where N is the number of initially introduced variables in phase space, we have 
D = 84— 56 — 2 x 6 = 16. This is in Hamiltonian formalism. In Lagrangian formalism, 
we have half of this number as dynamical variables. By dynamical variables, we mean 
those variables which obey differential equations containing accelerations. Taking a 
look on the complete action (|2.6| ) shows that the auxiliary variables A aX P are not 
within these variables. Therefore, Eq. fl3.20| ) says that after eliminating the redundant 



variables by using the constraints and gauge fixing conditions, we have 8 dynamical 



equations for g^s and £™'s. As we mentioned before, the total number of variables 
<7 M „ and is 4 times greater than the number of principle variable g^. Therefore, 
we conclude that the number of dynamical variables is 8/4=2 out of six g^ u . 

Notice that throughout our calculations concerning the constraint structure of 
the NMG model we have kept the cosmological constant term up to end. Hence, it 
is easy to find the constraint structure of the model without cosmological model just 
by putting A = 0. 

It should be noted that although the main structure of the constraints and the number 
of degrees of freedom is the same, the cosmological constant has a serious effect on 
the form of constrains [see eq.( |3.7[ )l as well as the final Hamiltonian, after imposing 
the derived forms of Lagrange multipliers in the total Hamiltonian ( |2.13| ). Therefore, 
it can be expected that the dynamics of the remaining physical degrees of freedom 
in the reduced phase space is affected deeply by the presence of the cosmological 
constant. Specially, particular solutions of the equations of motion appear when 
A ^ which are forbidden in the absence of cosmological constant. However, our 
purpose in this paper is not investigating the properties of particular solutions, which 
stand beyond studying the constraint structure of the system, although it could be 
interesting in turn. 

4. Concluding Remarks 

In this paper we studied the Hamiltonian structure of new massive gravity model. 
This is a complicated model in 3d gravity that contains higher order derivatives as 
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well as higher than quadratic terms. We used some form of the Palatini formalism 
in which combinations of Christoffel symbols, i.e., the variables are used as in- 
dependent variables while their relation with derivatives of the metric is imposed as 
additional conditions in the Lagrangian using the auxiliary variables A aX/3 . 
As is expected, the system is highly constrained with 36 primary and 26 secondary 
constraints, where 56 of them are second class and 6 are first class. This classi- 
fication makes the system suitable to be studied more carefully in the context of 
constrained systems. For example, one may be interested in finding the generator 
of gauge transformations in terms of first class constraints and study more carefully 
the gauge symmetries of the system. Moreover, if one decides to fix the gauges, 
one needs to know carefully which gauge fixing conditions should be imposed. As 
is well-known [16] and [18], for this proposes it is necessary to know the constraint 
structure of the system. 

We showed, finally in phase space there remain 16 physical variables. This means that 
there are 8 dynamical Lagrangian variables composed of the metric and Christoffel 
symbols. If one eliminates the Christoffel symbols, there remain 2 dynamical degrees 
of freedom out of 6 components of the metric. This conclusion is in agreement with 
the result that the NMG model constitute 2 gravitons under linearization of the 
equations of motion. 

Acknowledgments: the authors thank M. Alishahiha for his useful 
comments and bringing their attention to NMG model. 

Note added. -When completing our paper we observed reference [19] in which the 
total number of 2 degrees of freedom is derived in another approach. 
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